Ground state properties of a three-site Bose-Fermi ring with a small number of atoms 
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We investigate a three-site ring system with a small number of quantum degenerate bosons and 
fermions. By means of the exact diagonalization of the Bose-Fermi-Hubbard Hamiltonian, we show 
that the symmetry of the ground state configuration is a function of both the boson-boson and the 
inter-species interaction in the system. The phase diagram of the system, constructed by computing 
the exact two-body spatial correlations, reveals nontrivial insulating phases that exist even in the 
strong bosonic tunneling limit and for incommensurate filling of bosons. These insulating phases 
are due to the inter-species interactions in the system and are not necessarily accompanied by the 
suppression of the particle number fluctuations. 

PACS numbers: 



I. INTRODUCTION 

The ability to load ultracold quantum degenerate gases 
into periodic or quasi-periodic potentials of different ge- 
ometries, created by optical lattices [l| has opened the 
way to realization of well known condensed matter sys- 
tems with atoms of fermionic [3[ and bosonic character 
[H]. Such model systems allow us to gain deeper un- 
derstanding of the fundamental problems in many body 
physics beyond the weak coupling, mean-field regime, 
such as transition between superfluid and Mott-insulator 
In the strongly interacting regime, new phenomena 
can emerge, such as quantum magnetism Q and the su- 
persolid phase [3] ■ In general, these emergent phenomena 
belong to the class of quantum phase transitions and their 
mechanisms are strongly dependent both on the geome- 
try of the lattice and the interaction between the atoms 
loaded into the lattice. The theoretical and experimental 
studies of these effects in large many-sites lattice systems 
are important both from the fundamental and applied 
points of view. For example, the physics of ultracold lat- 
tice systems where frustration can occur, is fundamental 
for the development of fault tolerant systems and is im- 
portant for quantum computing Q. 

On the other hand, the possibility to address atoms 
in single sites 10] and to create mixtures of bosons 
and fermions |ll| draws the attention to the study of 
finite, inhomogeneous systems (l2j with complex inter- 
species interactions. The study of such two- or three- 
dimensional systems enables understanding of basic 
physics of more complicated systems, such as interaction- 
dependent properties of ground states (such as frustra- 
tions) and the role of symmetry breaking in transitions 
between different quantum phases. Small, strongly cor- 
related systems may also lend themselves to the process 
of controlled quantum state preparation and manipula- 
tion making it potentially useful in quantum informa- 
tion , q uantum-limited measurements, and atomtronics 

[Ski. 

In this paper we consider a minimal finite two- 
dimensional lattice model, namely a three-site ring of 
a Bose-Fermi ultracold mixture. Such a three-site sys- 



tem may be realized experimentally by engineering mag- 
netic microtraps on an atomic chip, or by combining a 
harmonic potential with a triangular or Kagome lattice, 
as suggested in With the small number of atoms, 

this system lends itself to the Bose-Fermi-Hubbard model 
solvable by means of direct diagonalization. We con- 
sider the ground state of the system and investigate how 
the admixture of fermions leads to various phases, de- 
pending on the filling factor and inter-species interaction 
strength. In particular, we consider unusual insulating 
phases resulting from the inter-species interaction, which 
is connected to the existence of macroscopic self-trapping 
states in the mean-field regime (l4| . 



II. THE MODEL AND THE GROUND STATE 
CONFIGURATION 

Based on the standard Bose-Fermi Hubbard model 
(see, e.g., 12]), the Hamiltonian for the three-site ring 
can be written as follows: 



H = H b + H f +H bb +rl bf , 



where. 
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Here w (b ) are the creation (annihilation) operators 
for the bosons and f' (f ) the creation (annihilation) 
operators for the fermions; the number operators are: 
= , £ e {b,f}- The nearest neighbour tunneling 
coefficient is T^, the intra- and inter-species interaction 
strengths are U bb > and U b f, respectively. 
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FIG. 1: (Color online) Schematics of the on-site state config- 
urations of three bosons (blue circles) and one fermion (red 
circle), with the lowest energy. The figures correspond to the 
states: (a) |1)+, (b) |2) + , (c) |3) + , and (f) |4)+ for repulsive 
inter-species interaction, and (d) (e) |2)_, (c) |3)_, and 

(f) |4)_ for inter-species attraction. 



In general, one can write a state vector of the system as 



follows: |*) = |*)/( 
where, n[ 



= N f , 



and n\ 



N h 



In the case of inter-species repulsion and Nf = 1, the 
states corresponding to the lowest energy levels are the 
ones with the following on-site particle distributions: 

|1)+ -> 11,0,0)/ ®|0, 3, 0) 6 , 

|2)+ -> |1,0,Q)/®|0,1,2) 6 , 

|3)+ -> 10,0,1)/® |0,2,1) 6 , 

|4)+ -> |1,0,0)/®|1,1,1) 6 , 

where |1) + , |2) + and |3)+ have six-fold degeneracy, while 
|4) + is a triplet. These states are schematically shown in 
Fig. [T] and the corresponding energies are plotted in Fig. 
1(a). 

For attraction and JV> = 1 (see Fig. [T]and Fig. [1(b)), 
the lowest-lying energy states are as follows: 

|1)_ -> |0,0, l>y (8) |0,0,3> 6 , 

|2)_ -+ |1,0,0)/®|2,0,1) 6 , 

|3)_ -»• |1,0,0)/®|1,0,2) 6 , 

|4)_ -> 10,1,0)/® |1, 1,1) 6 , 

where the states with symmetry |1)_ and |4)_ are triplets 
and |2)_ and |3)_ have six-fold degeneracy. Similar 
ground state configurations occur at incommensurate fill- 
ing, as shown in Fig. [3l 

As the bosonic interaction strength changes, the 
ground state of the system is changing symmetry. The 
hierarchy of the lowest lying energy levels correspond- 
ing to different states is shown in Fig. [2] for the com- 
mensurate and in Fig. 3] for incommensurate filling of 
bosons. In the case of commensurate filling, for a fixed 
value of the inter-species interaction and growing Ubb, the 
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FIG. 2: (Color online) The energy of the states with the 
largest contribution to the ground state as a function of Ubb 
for inter-species (a) repulsion and (b) attraction with the fixed 
magnitude \Ubf\ = 10. States with different symmetries are 
marked with the points. The unmarked states in (a) and (b) 
are |3)+ and |3)_, respectively. 



ground state structure evolves from a mixture of degen- 
erate states (a) and (b) to the state (/) (in Fig[T]) in the 
case of inter-species repulsion, and from (d) to (e) and (/) 
in the case of inter-species attraction. Similarly, in the 
case of incommensurate filling, the ground state evolves 
from a mixture of states (a) and (b) to (c) and to the 
mixture of degenerate states (c) and (/) (see Fig. [3]) for 
inter-species repulsion, and from (d) to (e), to the mix- 
ture of states (c) and (/) for the attractive inter-species 
interaction. 

The energy spectrum of this small scale system ex- 
hibits a gap, A, given by the energy difference between 
the ground state and first excited manifold [see Fig. [5] 
(a)]. For the commensurate filling of bosons this gap 
opens up at sufficiently strong boson repulsion for both 
repulsive and attractive inter-species interaction, as well 
as for Ufb = [see Fig. [5Jb)], and its magnitude is pro- 
portional to the interaction strength. It is reminiscent of 
the gap in the excitation spectrum of bosonic systems in 
lattices that indicates superfluid (SF) to Mott-insulator 
(MI) transition. Indeed, using the exact diagonaliza- 
tion of the Hamiltonian ([T]) and extracting characteristic 
behavior of tunneling correlations and particle number 
fluctuations, we can construct the phase diagram of the 
ground state and examine the transition to the insulating 
states in our small-scale Bose-Fermi system. 



III. PHASE DIAGRAM 

The typical characteristic quantity for different quan- 
tum phases of an ultracold quantum degenerate gas in 
a the lattice system is the spatial (tunneling) corre- 
lation for the bosonic component % and the tunnel- 
ing correlation of the fermionic component 77/ , given by 

% = l(£li + i>lM vg) £ G {6, /}, where n{ vg = N f /3. The 
boson tunneling correlation can be used to construct the 
phase diagram of the model, depending on the interac- 
tions. In a bosonic lattice this quantity tends to zero 
when the system is in the Mott-insulator regime, and ap- 
proaches one when the system is in the superfluid regime. 
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FIG. 3: (Color online) Schematics of the on-site state con- 
figurations of four bosons (blue circles) and one fermion (red 
circle), contributing to the lowest energy band. The panels 
(a), (b), (c), and (f) correspond to repulsive inter-species in- 
teraction and (d), (e), and (f) - to inter-species attraction. 




FIG. 4: (Color online) The energy of the states with the 
largest contribution to the ground state as a function of Ubb 
for inter-species (a) repulsion and (b) attraction with the fixed 
magnitude of inter-species interaction \Utf\ = 10. Energies of 
states with different symmetries shown in Fig. [3] are marked 
with the corresponding letters. 



Similarly, in our small, finite, mixed-species system, this 
quantity can be used to delineate between an insulating 
and a superfluid state. The fermion tunneling correla- 
tion measures the mobility of the fermion in the system. 
When its value is close to one, the fermion is mobile, and 
when it is zero, the fermion is pinned to one of the sites. 
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FIG. 5: (Color online) (a) Structure of the energy spectrum 
vs. Ubf for commensurate filling of bosons (Ubb = 1, iV& = 3). 
(b) The dependence of the gap, A = E\ — Eo on the intra- 
and inter-species interaction strengths for N b = 3. 



A. Commensurate boson filling 

As a start, we analyze the components of the ground 
state for the system with a commensurate number of 
bosons, Nf — 1 (n[ vg = 1/3), and for the fixed inter- 
species interaction strength \Ubf\ — 10. The boson tun- 
neling correlation for Nh = 3,9 is shown in Fig. [6] Due 
to the small number of particles, there exists a distinct 
crossover region between the SF (rjh = 1) and insulating 
(rjb = 0) regimes [l3j], where < r]b < 1. When this quan- 
tity is closer to either of the two extreme values, we will 
refer to it as a SF- or Mi-enhanced regime, respectively. 
It can be seen that the system exhibits a rich phase di- 
agram with an asymmetric behavior depending on the 
sign of the inter-species interaction and the number of 
bosons. 

For the case of attraction between bosons and 
fcrmions, Ubf < 0, a new insulating phase correspond- 
ing to the ground state |1)_ (Fig. [TJ d) appears. It 
is strikingly different from the regular Mi-like insulat- 
ing state in the pure bosonic system, |4)± (see Fig. [TJ), 
which appears at Ubf = and dominates the phase dia- 
gram for larger inter-species interaction strength. In this 
interaction-induced insulating phase the bosonic occupa- 
tion of the ring sites is strongly unbalanced, in analogy 
with the mean-field macroscopic self-trapped states of a 
Bose- Fermi mixture fl3| . As one increases the repulsion 
between bosons, a new SF enhancement region appears 
and then the ground state of the system is once again 
dominated by the regular insulating state |4)_. With 
the increasing inter-species interaction and the number 
of bosons, the SF enhancement region breaks into sev- 
eral filaments separated by the insulating states |2)_ that 
appear due to interaction between bosons and fermions. 
This change in the structure of ground state can also be 
followed in Fig. El (b). 

In the case of inter-species repulsion, Ubf > 0, the 
system is superfluid for small \Ubf\- In Fig. Ufa) one 
can clearly see that the ground state in this region is 
superposition of degenerate states |1)+ and |2) + , which 
is a typical sign of frustration. Insulating regions appear 
for larger inter-species interaction and are dominated by 
the |2)+ state [see Fig. [B](b)]. For large Ubb transition 
to the regular insulating state |4) + occurs. 

For both signs of inter-species interaction, the behav- 
ior in the regions where Ubb ~ Ubf is SF-enhanced. The 
effect of increasing the number of bosons is to scale up 
the regions of SF behaviour and introduce additional in- 
sulating regions due to inter-particle interactions. The 
appearance of the gap in the energy spectrum (Fig. [5]) 
correlates exactly with the insulating regions in the phase 
diagram. 

While the phase diagrams in Fig. (a,b) are based 
on the behaviour of the bosonic fraction, it is useful to 
examine the tunneling correlation of the fermion in the 
system. As seen in Fig. [S^c), in the region corresponding 
to the regular bosonic insulating state, (i.e. for \Ubb\ ^ 1) 
and for low inter-species coupling, the fermions are free 
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FIG. 6: (Color online) (a,b) Density plots of the boson tunnel- 
ing correlation, i] b (Ubf,Ubb), for (a) N b — 3 and (b) N b = 9. 
(c,d) Density plots of the (c) fermion tunneling correlation, 
t)f(Ubf,Ubb) and (d) boson number fluctuations, a(Ubf,Ubb) 
for Nh = 3. 



to hop between the ring sites, which can also be deduced 
from the symmetry of the |4)± state [see Fig. [Tj. In 
the mean-field picture [l4| this behavior reflects the fact 
that the effective interaction-induced potential seen by 
the fermion is weak and completely symmetric. As one 
increases |£4/|, the fermions localize and no tunneling is 
possible. In contrast to the regular insulating phase, the 
bosonic insulating phases that arise purely due to the 
inter-species interaction and correspond to symmetry- 
broken states (e.g., |1)_ or |2}+), naturally give rise to 
the regions of suppressed tunneling of the fermion. 

To characterize the interaction-induced insulating 
phases further, we look at the boson number fluctuations 
in the system given by a — \J (n\n\) — (h\) 2 /n 



avg ) 



see 



Fig|6fd). As expected from the MI behavior, the fluc- 
tuations in the regular insulating phase (with the |4)± 
symmetry) approach zero as one increases the strength 
of repulsion between bosons. In contrast, for attractive 
inter-species interaction the interaction-induced insulat- 
ing region at low values of Ubb is dominated by fluctu- 
ations, which can be taken as a signature of the new 
insulating state. 



B. Incommensurate boson filling 

It has been noted (see, e.g., [H, [Hj]) that for a few- 
particle pure bosonic system at incommensurate filling a 
small superfluid fraction is always present, therefore no 
insulating state can occur. We find that this is indeed 
the case in our system with no fermions (or Ubf = 0). 
The insulating phases in this case occur only in the pres- 
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FIG. 7: (Color online) Density plot of the (a) boson and (b) 
fermion tunneling correlation for Nf = l,Nb ~ 4. 
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FIG. 8: (Color online) Density plot of the (a) boson and (b) 
fermion tunneling correlation for Nf = 2,Nb = 10. 



ence of fermions, with non-zero inter-species interaction 
strength. For the repulsive inter-species interaction the 
insulating phase at large Ubb is suppressed, as shown in 
Fig. [Tfa), which is typical for Nb — 3m + 1 bosons with 
m being a positive integer. In contrary, the interaction- 
induced insulating region is very prominent for small Nb 
and is broken up by regions of enhanced superfluidity for 
larger Nb. For inter-species attraction both the regular 
MI- like insulating region at large \Ubb\ and the insulat- 
ing region at lower values of Ubb are retained due to the 
inter-species interaction. As in the case of commensurate 
filling, at large Ubb the new insulating domains that ap- 
pear due to the inter-species interaction are dominated 
by the ground state configuration (/) in Fig. [3J In gen- 
eral, both for incommensurate and commensurate filling 
of bosons the boundaries of the insulating regions are 
pushed to larger interaction strength values as one in- 
creases the number of bosons. 

In striking difference from the commensurate case, the 
tunneling of fermions is almost completely suppressed in 
the attractive inter-species interaction region, Ubf < 
[see Fig. EJb)]. For repulsive inter-species interaction, 
Ubf > 0, the fermion gains mobility as one increases the 
boson interaction strength, Ubb- This is due to the fact 
that, for larger Ubb and repulsive inter-species interac- 
tion, the insulating phase is dominated by the state (c) 
in Fig. [3J [see also Fig. Fig. @] (a)]. In this case the 
fermion is always able to hop between the sites with lower 
bosonic occupation numbers. In contrast, in the case 
of inter-species attraction the insulating phase is domi- 
nated by the state (/) in Fig. [3] [see also Fig. 2] (b)] and 
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the fermions are pinned to a site with the highest occu- 
pation of bosons. For weaker inter-species interaction, 
Ubf < Ubb-, the fermion is delocalized as in the commen- 
surate case. 



C. Role of the fermion filling factor 

So far, we have considered one fermion interacting with 
Nb bosons. Therefore, the results presented above are 
also applicable to a mixture of two bosonic species (see, 
e.g., [121) w hh a single atom in one of the components. 
In the case of Nf ^ 1, the Bose- Fermi- Hubbard model 
is still valid for the fermion filling factor less than 1 
(the system with no fermions is equivalent to the sys- 
tem with three fermions), and the fermion statistics in- 
fluences the ground state configuration. In particular, 
the system is particle-hole symmetric for 1/3 and 2/3 
filling of fermions. This fact is reflected in the behavior 
of the phase diagram, so that the case of 1 /3 fermion fill- 
ing with repulsive inter-species interactions corresponds 
to the case of 2/3 fermion filling with attractive inter- 
species interaction. For example, in the case of 2/3 filling 
of fermions the interaction- induced insulating regions ap- 
pear for repulsive rather then attractive interaction be- 
tween species, as seen in Fig. [H Other characteristic 



properties of the system, such as the behaviour of the 
fermion tunneling correlations and bosonic number fluc- 
tuations, are qualitatively the same as for Nf = 1. 

IV. CONCLUSIONS 

In conclusion, we have analyzed the ground state of 
a small-scale system of quantum degenerate bosons and 
fermions in a three-site ring configuration. We have re- 
stricted the consideration to the fermion filling factor less 
or equal than one, which as allowed us to employ a stan- 
dard Bose-Fermi-Hubbard Hamiltonian. By examining 
the tunneling correlations and particle fluctuations in the 
system, we have found that the system admits mobile and 
insulating states that are analogous to the superfluid and 
Mott-insulator states in infinite lattices. The novel insu- 
lating states identified in this small-scale system for both 
commensurate and incommensurate filling of bosons, are 
purely due to the inter-species interactions, and can be 
controlled by controlling the interaction strengths and 
the number of fermions injected into the system. 
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